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The POPOP4 code for converting secondary gamma-ray yield data to 
multigroup secondary gamma-ray production cross sections and the POPOP4 
library of secondary gamma-ray yield data are described. Recent results 
in the "testing" of uranium and iron data sets from the P0P0P4 library 
are given. The data sets were tested by comparing calculated secondary 
gamma-ray pulse-height spectra with spectra measured at the ORNL TSR-II 
reactor . 


One of the primary considerations in the design 
of shields for reactors on spacecraft is the dose 
rate due to secondary gamma rays produced by neutron 
interactions in the fuel and in the structural and 
shielding materials. The importance of secondary 
gamma rays is illustrated by the calculated dose 
distribution shown in figure 1 for a three-cycle W- 
LiH shield for a SNAP-8 ZrH reference reactor (ref. 

1) . The design objectives were to limit the dose 
rate on the outer surface of the shield to a maximum 
value of 1.55 x 10 3 mrem/hr and to minimize the 
weight of the spherical shield. Note from the figure 
that before the constraint dose is obtained, the 
secondary gamma-ray dose rate predominates. These 
dose rates were calculated by coupling the transport 
of the neutrons and secondary gamma rays in a single 
discrete ordinates calculation. One of the require- 
ments for a coupled calculation is accurate secon- 
dary gamma-ray production cross sections (SGRPXS f s). 
To satisfy this requirement, the shield designer 
must have a readily accessible library of secondary 
gamma-ray yield data for the elements used in shield 
design, the means to convert the yield data to 
SGRPXS 1 s in a format for use with the standard trans- 
port codes, and confidence in the accuracy of the 
data. 

P0P0P4, a FORTRAN-IV code, was written to con- 
vert secondary gamma-ray yield data as found in the 
literature for (n,y) , (n,xy) , etc., reactions to any 
required neutron-gamma multigroup energy structure 
(ref. 2). If the yield data are in terms of gamma- 
ray intensity per neutron induced reactions, the 
code multiplies the converted multigroup yields by 
input multigroup neutron-reaction cross sections to 
give the SGRPXS ? s. [The neutron reaction cross sec- 
tions are obtained from codes such as GAM-II or 
XSDRN (refs. 3 and 4).] P0P0P4 sums the SGRPXS's 
for the various neutron- induced reactions to give 


total SGRPXS's for the nuclide of interest. Using 
codes such as the Sample Simple Coupling Code 
(ASSCC) (ref. 5), the P0P0P4 cross sections are 
coupled with P^ neutron and gamma-ray cross-section 
sets for use in coupled neutron-gamma transport cal- 
culations. Discrete ordinates codes such as ANISN 
(ref. 6) and DOT (ref. 7) and Monte Carlo codes 
such as MORSE (ref. 8) use multigroup coupled cross 
sections produced as described above. 



Dose Distribution in Snap-8 Reactor with Three -Cycle W-UH Shield. 


FIGURE l.-Dose Distribution in SNAP-8 Reactor With 
Three-Cycle W-LiH Shield. 
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A compendium of neutron-induced secondary 
gamma-ray yield and cross-section data has been 
compiled for use with P0P0P4 • This collection of 
data is known as the P0P0P4 library (ref . 9) . In- 
cluded in the library are capture and inelastic- 
scattering yield data sets from the United Nuclear 
Corporation publications (refs. 10, 11, and 12), 
Maerker's and Muckenthaler *s measured secondary 
gamma-ray yields due to thermal-neutron captures in 
the elements found in soils, concretes, and struc- 
tural materials (ref. 13), the yields for thermal- 
neutron captures in a large number of elements as 
compiled in the Nuclear Data publications (refs. 14, 
15, and 16), the isotropic components of the gamma- 
ray production cross sections for Na, Mg, Cl, K, and 
Ca reported by Drake et al. (ref. 17), and many 
other (n,y) and (n,xy) data sets from the literature 
and from the private files of contributors. At pre- 
sent there are 223 data sets in the library for 79 
elements or nuclides. The library is available on 
magnetic tape from the Radiation Shielding Informa- 
tion Center (RSIC) . [P0P0P4, the Sample Simple 

Coupling Code, and P0P0P4 Library Tape Maker (ref. 

5) - a code to make or update the P0P0P4 library 
tape - are also available from RSIC.] 

One of the most important objectives of the 
P0P0P4 project is the "testing" of secondary gamma- 
ray yield data for use in shielding calculations. 

The data are tested by comparing measured secondary 
gamma-ray pulse-height spectra with spectra calcu- 
lated using P0P0P4 multigroup SGRPXS's prepared with 
data sets from the P0P0P4 library. 

A series of experiments was conducted by Muck- 
enthaler et al . of ORNL's Neutron Physics Division 
to measure secondary gamma-ray pulse-height spectra 
resulting from the exposure of various slabs of 
reactor core, structural, and shielding materials to 
collimated neutron beams (refs. 13 and 18). The 
slabs were exposed to bare, Cd-filtered, and 10 B- 
filtered neutron beams emanating from the ORNL TSR- 
II reactor. A simplified illustration of the exper- 
imental configuration is shown in figure 2. Differ- 
ential pulse-height spectra due to secondary gamma 
rays produced by neutron- induced reactions in the 
slabs were measured at a point 20 ft from the slabs. 
A borated polyethylene filter was placed between the 
NaI(T&) detector and the slab to reduce thermal-neu- 
tron effects in the detector. Some materials 


required an additional LiH filter in front of the 
borated polyethylene to thermalize scattered neut- 
rons. These experiments provided the standards for 
testing the data sets from the P0P0P4 library. 


0*m.-DW« M-I42T 



FIGURE 2.-TSR-II Experimental Configuration for 

Measuring Secondary Gamma-Ray- Pulse-Height Spectra. 


A flowchart illustrating the general procedure 
for calculating secondary gamma-ray pulse-height 
spectra is shown in figure 3. SGRPXS f s were pre- 
pared by using P0P0P4 to convert the yield data sets 
being tested to a 27 neutron - 60 gamma energy group 
structure and, if the data were given in terms of 
secondary gamma-ray yields, to combine the data with 
spectrally weighted, 27 group neutron reaction cross 
sections calculated with XSDRN or GAM-II. ASSCC was 
used to couple the SGRPXS's, P 3 neutron cross sec- 
tions, and P 3 MUG (ref. 19) gamma-ray cross sections 
for the element of interest. Using multigroup rep- 
resentations of the TSR-II beams and the coupled 
cross sections as input, the one-dimensional dis- 
crete ordinates transport code ANISN was used to cal- 
culate the secondary gamma-ray angular flux emanat- 
ing from the slabs. All ANISN calculations were made 
using a four-term expansion of Legendre polynomials 
and a 96th order of angular quadrature. 
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CALCULATED AND MEASURED 
SECONDARY y-RAY PULSE 
HEIGHT SPECTRA 


Colculotional Procedure for "Testing" Secondary Gamma-Ray Yield Data 

FIGURE 3.-Caleulational Procedure for 'Testing' 
Secondary Gamma-Ray Yield Data. 

The FORTRAN- IV code LINFOLU (ref. 5) was used 
to correct the ANISN angular flux for the geometric 
attenuation and for the attenuating effect of the 
detector neutron shield (s) - giving the intensity of 
secondary gamma rays striking the detector within 
ehch gamma group. One-fifth of the intensity within 

each group was assumed to be the intensity of five 
equally spaced discrete gamma energies within the 
group. Using Maerker's experimentally determined 
response functions for the Nal(TJt) detection system 
(ref. 13), the intensities of the discrete gammas 
were folded with LINFOLD to give the calculated 
spectra. Plots of the calculated and measured spec- 
tra were produced as the means for comparison. 

Initial efforts in the data testing program 
were concentrated on the 23 Al, Fe, Ni , Cu, and Pb 
(n,y) and (n,n'y) data sets from refs. 11 , 12 , 13, 

14, and 16. Pictorial comparisons of the calculated 
and measured pulse-height spectra resulting from the 
exposure of slabs of these materials to the TSR-II 
neutron beams are shown in ref. 20. The comparisons 
provide a means of evaluating the validity of the 
data sets for various incident neutron energy ranges. 
Recent efforts have been directed to the testing of 
the U, W, and Ta data sets in the library and to the 
testing of a newly acquired Fe (n,y) data set. To 
illustrate the testing procedure, recent results for 
U and Fe are described below. 


Mucken thaler et al. measured the spectra resul- 
ting from the exposure of a 10.44.x 5.75 x 0.0304- 
in. -thick depleted uranium foil to the collimated 
TSR-II beams. Background was determined for each 
beam by measuring the spectra without the foil in 
place. The 235 U and 238 U number densities in the 
sample were 8.305 x 10 ~ 5 and 5.02 x 10 “ 2 nuclei/ 
barn cm, respectively. The foil was oriented as 
shown in figure 2 above. Since gamma rays are pro- 
duced by the fission capture, non-fission capture, 
and the inelastic-scattering reactions in both 23 5 U 
and 238 U, and since the SGRPXS 's depend on the res- 
onance characteristics of the nuclides, the spectra 
comparisons for the foil provided an interesting 
challenge. XSDRN was used to calculate 235 U and 
238 U P 3 neutron cross sections for the 27 neutron 
group structure. This calculation included the 
resonance self-shielding effects in both nuclides. 
Booth has shown that the secondary gamma-ray spectrum 
resulting from exposure of the foil to the Cd-filtered 
neutron beam minus 1.377 times the 10 B-filtered 
spectrum eliminates the effect of high-energy (fis- 
sion) neutrons (ref. 21). Sample comparisons for 
the (Cd-1,377 10 B) and bare neutron beams are shown 
in figure 4. The data sets used in the preparation 
of the SGRPXS 's for the transport calculations are 
listed in Table I. Calculation B is a bare beam 

calculation which included SGRPXS ’ s for all 235 U and 
238 U’ neutron induced reactions except the 238 U (n,y) 
reaction. SGRPXS's for all 235 U and 238 U neutron 
induced reactions were used in calculation D. The 
238 U SGRPXS 1 s for calculation D were prepared from 
data sets 928112, 928301, 928901, 925801, and 925804 
as indicated in Table I. Calculation A is a (Cd- 
1.377 1 °B) calculation which included SGRPXS's for 
all reactions except the 238 U (n,y) reaction, where- 
as, calculation C included SGRPXS T s for all 235 U and 
238 U reactions. A report describing similar tests 
of twelve uranium data sets from the P0P0P4 library 
is in the final stages of preparation (ref. 25) . 

Measured spectra resulting from the exposure 
of an approximately 5-ft-square by 1/16-in. -thick 
Fe slab to the TSR-II beams were used as the stan- 
dards for testing Fe data sets. Since the Cd-filter 
has a "cutoff" at 0.5 eV and a constant attenuation 
of 0.986 to neutrons above 1 MeV, and since the bor- 
on filter "cuts off" at approximately 10 keV and has 
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SECONDARY GAMMA-* RAY SPECTRA (counts/MeV/min * kW) 
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FIGURE 4. -Comparisons of Calculated and Measured 
Secondary Gamma-Ray Pulse-Height Spectra for a 
Depleted Uranium Foil. 

*Ordinate caption should read Secondary Gamma-Ray 
Pulse-Height Spectra (counts/MeV/min.kW) . 
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Table X. Uranium Secondary Gamma-Ray Yield Data 


Data Set 
Identification 
Number^ 

Uranium 

Isotope/ 

Reaction 

Source of 
Yield Data 

Remarks 

925101 

5 / (n , y) 

Ref. 12, p. 45 

Yields given for 11 neutron groups (3.7 x 10 2 £ E n £ 18.2 x 
10 6 eV) and 13 discrete gamma energies (0.0894 £ Ey £ 6.42 
MeV). Data used for the 235 U(n,Y) reaction in all calcu- 
lated spectra shown in Fig. 4. 

925301 

5/ (n ,n f ) 

Ref. 12, p. 52 

Yields given for 80 neutron groups (1.27 x 10 4 £ E n £ 18.017 
x 10 6 eV) and 11 gamma groups (0.0 £ Ey £ 6.50 MeV). Data 
used for the 2S5 U(n,n T ) and (n,2n f ) reactions in all calcu- 
lated spectra shown in Fig. 4, Assumed applicable to U 

in Calculations A, B, and C. 

925801 

5/(n,Y) F * 
and (n,2n') 

Ref. 22 

Thermal— neutron prompt and delayed fission capture yields 
given for 20 gamma groups (0.0 £ Ey ^6.50 MeV). These data 
assumed applicable to both 23 U ana U and to all incident 

neutron energies in Calculations A, B, and C. Data used in 
Calculation D for all neutron energies in 235 U SGRPXS and 
for En < 1.1 MeV in 238 U SGRPXS. 

925804 

5/(n,Y) p 

Ref. 22 

Delayed part of 925801 data (t > 1 sec). Data used only in 
Calculation D for En >1.1 MeV in 23 8 U SGRPXS. 

928113 

8/(n,Y) NF ! 

Ref. 23, 
pp. 51-65 

Yields given for 21 neutron groups (5.6 £ En £ 1.019 x 10 5 
eV) and 16 gamma groups (0.903 £ Ey £ 4.807 MeV). Data used 
in Calculation C for 230 U(n,Y) SGRPXS in applicable neutron 
energy range. 

928112 

B/(n,Y) NF 

Ref. 26 

Yields given for 37 neutron groups (0.0 £ En £ 11.08 MeV) 
and 117 gamma groups (0.0 £ Ey £ 5.85 MeV). Data used in 
Calculation C for E n < 5.6 eV and E n > 1.019 x 10 5 eV. 
Data used in Calculation D for En < 1.1 MeV. 

928301 

8/ (n,n f ) 
and (n,2n*) 

Ref. 12, p. 54 

Yields given for 79 neutron groups (0.337 £ E n £ 18.02 MeV) 
and 11 gamma groups (0.25 £ Ey £ 6.5 MeV). Data used in 
Calculation D for En < 1.1 MeV. 

928901 

8/(n,Y) NF , 

(n,Y)p» and 

(n ,n* ) 

Ref. 24, p. 19 

SGRPXS given for 15 neutron groups (0.485 £ En £ 15.3 MeV) 
and 12 gamma groups (0.50 £ Ey _< 6.5 MeV) . X— sects used m 

Calculation D for 238 U(n,Y)flF» ( n >Y)prompt F> and ( n » n? ) 
reactions in applicable neutron energy range. 


(n,Y) NF is the nonfission capture reaction and (n,Y)p is the fission capture reaction. 

Each data set in the P0P0P4 library is identified by a unique six-digit identification number 

Table 11. Iron Secondary Gamma-Ray Yield Data 


Data Set 
Identification 
Number 



Source of 
Yield Data 

Ref. 13, 
p . 67 


Remarks 

Thermal -neutron capture yields 
given for 80 discrete gamma 
energies (1.24 £ Ey £ 9.30 MeV) . 
These data assumed applicable to 
all incident neutron energies in 
calculation labeled n ORNL-4382 


(n,n') 


Ref. 26 Yields given for 40 neutron 

groups (0.0 £ E n <_ 1.0026 MeV) 
and 200 gamma groups (0.0 £ Ey 
£ 10.0 MeV). Data used in cal- 
culation labeled M JKY, JEW . ..." 

Ref. 12, Yields given for 53 neutron 

p. 35 groups (0.84 £ E n £ 18.1 MeV) 

and 10 gamma groups (0.50 £ Ey 
£10.0 MeV). Data used in both 
calculations shown in Fig. 5. 



a constant attenuation of 0.75 to neutrons above 1 
MeV, the spectra from the (Cd- 10 B) neutron beam pro- 
vide a test of the epithermal Fe capture yield data 
(ref. 18). Two calculated spectra resulting from 
exposure of the slab to the (Cd- 10 B) beam are com- 
pared with the measured spectrum in figure 5. The 
data sets used in the preparation of SGRPXS for the 
calculations are identified in Table II. Neutron 
cross sections for the calculations were calculated 
with GAM-II. 


Results of these studies have proven that the 
conversion technique used in P0P0P4 is satisfactory 
for the conversion of yield data to a required mul- 
tigroup energy structure. Since the SGRPXS' s are 
only as accurate as the yield data and the neutron 
reaction cross sections, the neutron cross sections 
must be properly weighted prior to being used. The 
growth envisioned for the P0P0P4 library depends on 
the continued support of individuals who are willing 
to share their efforts with others. 
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